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Upcoming imaging surveys such as the Large Synoptic Survey Telescope will repeatedly scan 
large areas of sky and have the potential to yield million-supernova catalogs. Type la supernovae 
are excellent standard candles and will provide distance measures that suffice to detect mean pair- 
wise velocities of their host galaxies. We show that when combining these distance measures with 
photometric redshifts for either the supernovae or their host galaxies, the mean pairwise velocities of 
the host galaxies will provide a dark energy probe which is competitive with other widely discussed 
methods. Adding information from this test to type la supernova photometric luminosity distances 
from the same experiment, plus the cosmic microwave background power spectrum from the Planck 
satellite, improves the Dark Energy Task Force Figure of Merit by a factor of 1.8. Pairwise velocity 
measurements require no additional observational effort beyond that required to perform the tra- 
ditional supernova luminosity distance test, but may provide complementary constraints on dark 
energy parameters and the nature of gravity. Incorporating additional spectroscopic redshift follow- 
up observations could provide important dark energy constraints from pairwise velocities alone. 
Mean pairwise velocities are much less sensitive to systematic redshift errors than the luminosity 
distance test or weak lensing techniques, and also are only mildly affected by systematic evolution 
of supernova luminosity. 
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I. INTRODUCTION 

Measurements of the distance-redshift relation of type la supernovae (SNela) firmly established contemporary 
accelerated cosmological expansion [lj, |2[ and SNIa distances remain one of the most promising probes of dark energy. 
To determine whether the accelerated cosmological expansion is caused by an ubiquitous dark energy or large-scale 
deviations from general relativity, it is necessary to measure both the expansion of the universe and the dynamics of 
the structure formation [3l-[l4j]. The SNIa luminosity distance test provides information about the expansion rate of 
the universe, but does not provide information on structure formation (though SNIa magnifications may achieve this 
in the future, see Refs. |15l-lr7|). Peculiar velocities are related to densities through a continuity equation, so peculiar 
velocity statistics provide one avenue to study the growth of cosmic structure (e.g., [II]). The most well-explored 
option for probing the peculiar velocity field is via redshift-space distortions imprinted on the galaxy power spectrum 
(e.g., j!9l42l| ). Peculiar velocities may be detectable with future microwave experiments via the kinetic Sunyaev- 
Zeldovich effect [22T - [24| and from large samples of SNela with spectroscopic redshifts 25]. In this paper, we examine 
the possibility of utilizing the mean pairwise velocity statistic, measured from SNela in a large photometric survey, 
to constrain dark energy. 

Two well-studied statistics derivable from a sample of line-of-sight peculiar velocities are the velocity correlation 
function and the mean pairwise velocity [26l . [27j . The former is a two-point statistic expressing correlations in the 
peculiar velocities of objects are as a function of their separation. The mean pairwise velocity is a measure of the typical 
relative velocity of objects at a given separation. Peculiar velocities are sensitive to both the rate of structure growth 
in the universe and the rate of expansion of the universe. Therefore, peculiar velocity measurements on cosmological 
scales may constrain the dark energy that drives cosmological acceleration and quenches late-time structure growth. 

Traditionally, the bulk flow velocity has been measured by coupling measured galaxy redshifts with local distance 
indicators such as the fundamental plane of early -type g alaxies (28l . |29| , the Tully-Fisher relation [3fj| - l32l | , or surface 
brightness fluctuations [33|. More recent studies |34l - [3a | have measured significant bulk flows on scales of 100 Mpc. 
Radial velocity measurements have been used to reconstruct the velocity and density fields [37j . Reconstruction 
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methods provide a way to test the gravitational instability theory and to measure the bias between the galaxy 
and mass density fields. Such studies are limited to the relatively local Hubble flow (z < 0.1), primarily because 
a constant fractional error in distance corresponds to a larger velocity interval at higher redshifts, an error that 
eventually overcomes the signal. 

Type la supernovae, in contrast, are well-calibrated standard candles, and at cosmological distances SNela are 
more reliable distance indicators than those previously used for measuring peculiar velocities. Indeed, the dipole and 
quadrupole moments of the local bulk flow velocity have been measured to higher precision with the current data set 
of a few hundred SNela than with reconstructions based upon catalogs of many thousands of galaxies [25J . 



SNela that are physically near each other exhibit coherent motion as they are influenced by correlated density 
structures. Therefore, the errors in the luminosity distance measurements of pairs of SNela should be correlated at 
low redshift (z ^ 0.1). Ignoring this correlation can lead to systematic biases in the determination of dark energy 
parameters [38-40]. Alternatively, one can treat these correlated shifts in luminosity distance as "signal," because 
peculiar velocities depend upon cosmological parameters. This signal has led to useful, independent constraints 
on the low-redshift normalization of the matter power spectrum, as, and the total matter density, VL m [il], [42J. 
Unfortunately, direct measurements of the velocity correlation remain limited to relatively low redshifts. Even in an 
optimistic scenario of measurements of one million SNela, all with full spectroscopic follow-up, the velocity correlation 
can only be measured to a redshift of z ~ 0.5 [Tl| . 

In contrast to peculiar velocity correlation measurements, mean pairwise velocity is a linear statistic so its errors 
vary more mildly with redshift. In this study, we show that it will be possible to obtain interesting cosmological 
information from mean pairwise velocities to a redshift of z — 0.9 in a large photometric survey of SNela, such as that 
planned for the Large Synoptic Survey Telescope (LSST) which is anticipated to increase dramatically our current 
catalog of SNela, by a factor of nearly 1000. We demonstrate that such a measurement can provide dark energy 
constraints that complement luminosity distance measurements under optimistic, but reasonable, assumptions. The 
constraints from mean pairwise velocities are also useful because they may be estimated with relatively little additional 
observational effort beyond that already required to use SNela to map luminosity distance or to detect cosmic lensing 
magnification We find that combining the mean pairwise velocity measurements with distance measurements of 
SNela will sharpen constraints on the dark energy parameters compared to those inferred from luminosity distances 
alone. In particular, mean pairwise velocity constraints can improve the dark energy Figure of Merit from SNela 
as defined by the Dark Energy Task Force [43[ (DETF) by a factor of 1.8. We additionally demonstrate that mean 
pairwise velocities, being a differential statistic, are potentially much less sensitive to systematic errors than other 
commonly considered observational techniques. Ultimately, this property may make mean pairwise velocities one of 
the most practically useful probes of dark energy. 

Following the DETF, we describe the dark energy in terms of three phenomenological parameters: its current 
energy density tt\ and two parameters describing the redshift evolution of its equation of state, wq and w a , such that 
w(a) = wq + (1 — a)w a - The additional cosmological parameters upon which the velocity field depends are the large- 
scale normalization of the matter power spectrum Af, the power-law index of the primordial power spectrum ns, the 
Hubble parameter h, the curvature of the universe fifc, and the present-day matter density f2 m . In addition, we treat 
the photometric redshift (photo-z) dispersion, a z , as a free parameter with priors. We label our set of parameters p. 
We consider a fiducial cosmological model similar to the WMAP 5-year results [44]: Af = 2.0 x 10~ 9 , ns — 0.95, 
h = 0.71, fl k = 0, fl m = 0.25, fl A = 0.75, w = -l, and w a = 0. 

The paper is organized as follows. In Section II, we describe our assumed SNIa survey specifications and review the 
estimation of supernova line-of-sight peculiar velocities from observed supernova brightnesses and redshifts. Section III 
describes a halo model calculation of the mean pairwise velocity as a function of cosmological parameters. Sections IV 
and V quantify various sources of systematic and statistical errors that impact SNIa pairwise velocity measurements 
respectively. We present our results for dark energy parameter constraints in Section VI, using two different sets of 
prior constraints. We also derive limits on systematic effects that must be obtained to have the resulting parameter 
bias be smaller than the calculated statistical errors. In Section VII, we summarize the kinds of observational efforts 
required to meet the prospects outlined in this paper, along with a brief discussion of the systematic error properties 
of mean pairwise velocities compared to other dark energy probes. 



Forthcoming large-scale imaging surveys such as LSST or the Panoramic Survey Telescope and Rapid Response 
System (PanSTARRS) [45M47I] will discover 10 4 to 10 6 SNela. These SNela may be observed with broadband pho- 
tometry with exposures spaced several days apart. To infer cosmological parameters from peculiar velocities, reliable 
distance measurements are needed. These will likely be obtained from a well-characterized subset of the supernovae 
discovered by any survey. The particular characteristics of this subset depend upon survey strategy and are difficult 
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to anticipate. 

For ease of comparison with published studies, we adopt survey specifications similar to what may be achieved 
with a survey similar to LSST. We assume a total of 3 x 10 5 SNela out to z = f.2, collected over a dedicated 
supernova survey region of 300 square degrees; this corresponds to a SNIa surface density of 1000 deg~ 2 . This 
number density corresponds to the "d2k" survey described in [48[ and such a dedicated survey may be undertaken 
as part of the science goals of the LSST [47]. We assume redshifts estimated using broadband photometry with a 
redshift-dependent, normally-distributed error of a z = <j z q(1 + z). DETF specifies an error range of <j z q = 0.01 for 
an optimistic scenario to a z o = 0.05 for a pessimistic scenario; in our parameter forecasts, we allow a z o to vary along 
with the cosmological parameters. Following Zhan et al. [lij], we model the SNela redshift distribution as 

d 3 n fexp(3.12z 21 ) - 1, z < 0.5, 

dfldzdt * j(exp(3.12z 21 ) - l)exp(-12.2(z-0.5) 2 ), z > 0.5. ^' 

To the extent that SNela are standardizable candles, photometric observations will yield a distance modulus [i and 
a luminosity distance dh via the usual relation 



^= 2 - 171n U^J +25 - (2) 

The luminosity distance is obtained from the cosmological redshift z via the definition 

f z dz' 

d L (z) = (1 + z)d c (z) = (1 + z)c J (3) 

where dc{z) is the comoving line-of-sight distance to a galaxy at redshift z, H(z) is the Hubble parameter as a function 
of redshift, and a geometrically flat universe has been assumed in the second equality. The evolution of the Hubble 
parameter, and thus the luminosity distance, depends on the assumed cosmological model. For a given supernova, 
its measured redshift is the difference between its cosmological redshift and the additional Doppler shift due to its 
line-of-sight velocity, 

Zmeas = z{jl) - — (1 + z(fi)), (4) 
c 

where its cosmological redshift z(p) can be obtained from its observed luminosity by inverting Eqs. ([2} and (|3|). The 
factor of (1 + z) in Eq. (j4j accounts for the cosmological redshift between the rest frame and the observation frame. 
For a given supernova with observed redshift and luminosity, its line-of-sight velocity can be obtained by rearranging 
Eq. dD into 

Czili) Corneas /_\ 

"los = — , (5) 

-L i ^rncas 

where we have replaced by z mcas in the denominator, which will always be a good approximation for objects at 
cosmological distances where the first term in Eq. (Q} is large compared to the second term. 

Traditional peculiar velocity estimates using other standard candles at cosmological distances have been hampered 
by errors in distance estimates, which propagate into errors in z(fx). For a galaxy with cosmological redshift z = 0.03, 
a 10% error in distance corresponds to an error in inferred cosmological redshift equivalent to a peculiar velocity of 
1000 km/s, with the size of the error increasing proportional to redshift for z < 1. Large-area supernova surveys offer 
two main advantages. First, supernovae are bright enough and good enough standard candles to provide convenient 
distance estimators out to z = 1 and beyond. Second, the anticipated large number of supernovae hold the promise of 
determining average distances far more precisely than individual distances, allowing precise determination of average 
velocity statistics from large catalogs of supernovae. Of course, realizing this promise requires controlling systematic 
errors in both distance and redshift observations to a high level, so that averages over large ensembles of SNela reflect 
the actual velocity statistic. Both systematic and statistical errors will be considered following the next Section, which 
outlines the application of the mean pairwise velocity statistic to supernova surveys. 



III. MEAN PAIRWISE PECULIAR VELOCITY 



The mean pairwise velocity v{r,a) at a comoving separation r and scale factor a = 1/(1 + z) is the average over 
all pairs at a fixed comoving separation of the relative peculiar velocity of the two galaxies projected along the line 



4 



joining them. That is, 

where is the peculiar velocity of supernova i and f is the unit vector in the direction of the separation of the 
two objects. The sum is over N(r) pairs at a given comoving separation r. (Note that the quantity which we write 
throughout this paper as u v(r, a)" is commonly written in the literature as "uy(r, a)" or ll V\%{r, a)." We use this 
notation to avoid potential confusion with subscript labels for individual galaxies that we use below.) 

The mean pairwise velocity for dark matter particles may be derived using the pair conservation equation [4^ |. 
However for galaxies, the pair conservation equation needs to be modified to account for evolution }27| . The resulting 
mean pairwise velocity for SNIa host galaxies with a comoving separation r at a mean scale factor a (assuming that 
the redshift difference between the two galaxies corresponds to a scale factor difference much smaller than a) can be 
written 

v(r,a) = --H(a)a^—b gal (a) 1 + ^ l{ray (7) 



where 



D 2 

fW ro ) = — £-/ dkksm(kr)P(k) (8) 







is the dark matter two-point correlation function, P(fc) is the dark matter power spectrum at wavenumber k, H(a) 
is the Hubble parameter at a given redshift, and D a is the linear growth factor as a function of time, normalized so 
that D a = 1 at z = 0. We also define the dark matter correlation function averaged over separations less than r to be 



e dm (r,a) = -| / dr' r' 2 e m {r\a). (9) 



We are interested in the large-scale limit, so we model the correlation function of supernova host galaxies using a 
deterministic linear bias relative to the dark matter b ga \{z), defined by 

^\r,a) = bl 1 (z)^(r,a). (10) 

The bias b ga \(z) in general varies with the galaxy separation, a variety of galaxy properties, and redshift |50j | . In 
the large-scale limit, scale-independent bias is a fairly good assumption. Following Ref. [48|, we model b ga i(z) as 
frgai(z) = 1-0 + 0.6z to obtain the fiducial value of the galaxy bias as a function of redshift. 

With future photometric surveys potentially detecting more than a billion galaxies, we can expect that the correla- 
tion function of samples of galaxies matching the SN hosts can be measured to percent-level accuracy or better. Thus 
the uncertainty in b ga \ will primarily be due to uncertainty in the cosmological parameters affecting the dark matter 
correlation function. We express b ga \ in terms of the galaxy and predicted dark matter correlation functions, and use 
this bias value in Eq. (J7J. 

Only the line-of-sight component of the velocity can be obtained from observations, while the mean pairwise velocity 
involves all three directional components of the velocity. We use the estimator for the mean pairwise velocity given a 
data set of line-of-sight velocities developed in Ref. [26| . Consider two galaxies i and j at comoving positions and 
Yj moving with peculiar velocities Vj and Vj . The radial component of velocities can be written as v\ = r j • Vj and 
Vj = Tj-Vj. Then an estimate for the pairwise velocity of the two galaxies i;??* is defined by (uf — up = v°j t r-(r i +ij)/2, 
where f is the unit vector along the line joining the two galaxies. If we now consider a catalog of line-of-sight galaxy 
velocities, minimizing % 2 between the actual pairwise velocities and the estimate of the pairwise velocity at a given 
separation r gives an estimator for the pairwise velocity Eq. (JTJ) based on the catalog, 

v^(r,a) = Ep ^~f )Pl \ (11) 
/—/pairs Pij 

where the sums are over all pairs i ^ j of galaxies at comoving separation r and — f • (?j + Vj)/2. Note that this 
form for the projection tensor is applicable in the flat-sky limit and breaks down for large angular separations; in 
particular it is zero if the two galaxies are in opposite sky directions. In this paper, we consider a model supernova 
survey of 300 square degrees in a compact sky region, and the expression given here is always valid. To extend 
the results here to a full-sky survey, or to survey patches which are separated by large angles, a more complicated 
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projection tensor must be used. The derivation is not conceptually difficult, but this will be deferred to future work 
giving more detailed estimates of signal-to-noise ratios for particular observing strategies. 

Equation is a function of the separation r between the two galaxies. To measure this distance, we must use 
the estimated locations of each galaxy; this is subject to errors which will be quantified in the next Section. The 
separation that is measured directly is the angle between two galaxies on the sky. This angle can be converted to the 
transverse component of the distance between the two galaxies using the angular diameter distances corresponding 
to their redshifts. 

The expression in Eq. (jlip is a very simple estimator which weights all pairs of velocities uniformly. A more careful 
analysis of real data would use, for example, a signal-to-noise weighting in the sum. This is not a major correction to 
the analysis in this paper, as we limit the sums in Eq. (fTTj) to pairs with separations smaller than 100 Mpc; at larger 
separations the signal becomes small. In principle, a signal-to-noise weighting can squeeze more information out of the 
data, using pairs with larger separations, but it does not qualitatively change our results. Our estimator is accurate, 
as we have shown explicitly in Fig. 4 of Ref. [23j . but suboptimal; an optimal estimator will somewhat improve the 
constraining capability of a velocity survey compared to the analysis here, so our estimator is conservative. 

As we discuss further in Section IIVI and Section IVI1 we can mitigate the influence of systematic redshift errors 
by considering a related projected statistic, where the mean pairwise velocity is taken as a function of the angular 
separation of the two galaxies rather than as a function of their three-dimensional separation. This is given by 

" dir t P(ir t \8,a)v(r,a), (12) 

where the line-of-sight comoving separation ir t = ^0(^2) — ^c(ai) and P(ir t \9,a) is the probability that a pair has 
line-of-sight separation 7r 4 given that it has an angular separation on the sky 9. We can write the three-dimensional 
separation r in terms of the angular separation 9 and 7r t as 



>d M {ay+irl (13) 



where 



cHq 1 ^ 1/2 sinh 



n k >o 



[n]! 2 d c {a)/{cH^ 

dM{a)= \d c (a), n k =0 ( 14 ) 

.ci^lfifel-^sin [\^ k \ x / 2 dc{a)/{cH^)) , n k < 

is the transverse comoving distance to scale factor a; here il k is the effective curvature density. f2fc = 1 — Cl m — CIa (see 
Ref. [5II ] for a lucid discussion of various distance measures in cosmology) . If the redshift difference is small compared 
to unity, ■K t H{zi) « c(zi — z\), though we compute the separation in full for all pairs. For a spatially flat universe, 
dhiio) — dc(a)- In our case, we always consider separations with r <§; cH^ 1 since the signal is only significant on 
these scales. We therefore always have dM{a) ~ dc(a) to good accuracy, and for simplicity we make this assumption 
throughout the rest of the paper and use comoving distances entirely. We consider pairs of galaxies with line-of-sight 
comoving separations up to a maximum value 7r max (in practice, we will measure redshift-space rather than comoving 
separations; we consider the impact of this in Section [IV A[) . The probability of a pair having line-of-sight separation 
7Tt given that it has an angular separation on the sky 9 is 

P(7r t |g,a)= 1 + f^ (r ;j (15) 
Jo dir t [l+^ l {r,a)} 

for 7r t < 7r max and P(n t \9,a) = for n t > 7r max . 

An estimator for v(9,a) from line-of-sight velocity data is easily obtained by substituting v est (r, a) for v(r,a) in 
Eq. (|12|). To compare with data, we bin this statistic in angular separation and redshift, putting each pair in the 
redshift bin corresponding to the mean photometric redshift of the two galaxies in the pair. In this manner all pairs 
are included regardless of binning; we have verified that our results remain similar when modest changes are made 
to projection and binning schemes. Note a correction for scatter in measured redshifts must also be included, as 
discussed below in the following Section. 

Changing the maximum separation 7r max considered in Eq. (I12[) will modify the signal-to-noise ratio in measuring 
the projected pairwise velocity. A larger 7r max increases the total number of pairs considered, but the signal-to-noise 
ratio for each pair decreases at larger ir t (as measurement errors remain approximately unchanged but signal strength 
decreases), so their contribution is small. For the purposes of this paper, we adopt a cutoff of 7r max = 100 Mpc, which 
captures the great majority of the pairwise velocity signal. As a test of this effect, we find that including pairs out to 
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separations two times larger only changes the signal-to-noise in measuring the projected pairwise velocity by around 
10%. Based on this, we conclude that including data from pairs with separations larger than 100 Mpc should give 
only minimal improvements in parameter constraints compared to those presented here. We also impose a minimum 
separation of 20 Mpc on the pairs we consider, to eliminate any systematic errors related to nonlinear effects. The 
mean pairwise velocity is a declining function as separation increases from 20 Mpc to 100 Mpc, as shown in Fig. [T] at 
smaller scales, it turns over and decreases in linear theory. 



IV. SYSTEMATIC ERRORS 



A. Photometric Redshift Errors 



Large imaging surveys will detect so many galaxies that it will not be feasible to obtain spectroscopic redshifts 
for the vast majority. We must settle for photometric redshift estimates determined from the fluxes measured in 
the various observed bands. These photometric redshifts will be less accurate than spectroscopic redshifts, and may 
have complex error distributions. Here we consider a measured redshift distribution described by a Gaussian of 
standard deviation a z centered at the true redshift of each object. We neglect a possible photometric redshift bias 
for two reasons: First, in realistic surveys this bias can be calibrated by comparison with a manageable number of 
spectroscopic SNIa observations [lT], EH, Sll- We emphasize that we utilize a normal distribution for definiteness, but 
a well-calibrated error distribution is what is necessary to proceed; errors need not be Gaussian in practice. Second, 
the expected level of photometric redshift bias is likely to be a small effect [H, H3l compared to the systematic errors 
in estimating distances that we consider below. As a result, we do not explicitly carry a bias through in the equations 
below, but we will present a test of the impact of a bias in photometric redshifts in Sec. IVII 

In contrast, the photo-z dispersion, <r Zl essentially smooths the estimated velocity distribution of the observed 
sample and propagates scatter into galaxy pair separations. The latter effect can cause not only a scatter in inferred 
cosmological parameter values, but also a systematic shift, which we calculate here. 

The mean pairwise velocity w(r, a), given in Eq. ([7]l. assumes that the three-dimensional separation r between 
the SNela or their host galaxy pairs are known accurately; however, there will be non-negligible errors in observed 
redshifts. Our simple normal-error model for the distribution of the photometric redshift z p , given a true redshift z, 
is 

P(z p \z,a z ) = -=L=exp[-(z-z p ) 2 /(2a 2 2 )]. (16) 

V 27T(J z 

We take the photo-z dispersion to be a z = a z a(l + z) with ct z0 ranging from 0.01 to 0.05 0, Hj| [H, El, Eil • We explore 
the sensitivity of our results to prior knowledge of a z0 in § IVII 

Using Eq. (|T6| and the expression H(z p )ir t = c(z P 2 — z p \) for the local Hubble expansion about each SNIa, where 
z p 2 — z p i is the photometric redshift difference between a pair of supernovae, we write the probability of obtaining 
the observed line-of-sight separation 7r b s for a given, true comoving line-of-sight separation 7Tt as 

P(*ob,\n,<r*) = J— y exp [-(7r obs - TT t ) 2 /{2aD] , (17) 
V 27T(J l 

where cv = \/2c<j z /H(z) [Hjjj . We assume that the photometric redshifts z p , although they include the effects of 
peculiar motions, give a better measurement of the galaxy line-of-sight separation than the cosmological redshifts 
z(/i), which must be determined via a distance measurement with uncertainties on the order of 10%; hence the line- 
of-sight positions of SNela are estimated using z p . The factor \/2 in relating ov to a z accounts for uncertainties in the 
positions of the two galaxies in a pair, which are added in quadrature. Combining Eqs. ([7]), ([121 . and (|17l) . we get an 
expression for the projected, mean pairwise velocity accounting for a significant dispersion in photometric redshifts, 

/•7T max P OO 

v(6,a\a^(a)) = / dirt / dTT ohs P(Tr t \9,a)P(TT ohs \Tr t ,a^(a))v((9 d c (a) 2 +n 2 y/ 2 ,a). (18) 
Jo Jo 

We propose using this statistic as a cosmological probe. We consider only positive values of 7r t , so we count each pair 
only once. This remains true if in some cases (due to errors) 77* scatters below zero (in which case the separation is 
positive when the two members of the pair are exchanged). 

Both Eq. (|17[) and the expression for a v are valid only when \z p2 — z p i \ <C 1; however, we should always be in this 
limit. The maximum true separation we consider, 7r max = 100 Mpc, corresponds to a redshift difference ranging from 
0.024 to 0.042 as z ranges from to 1; photo-z errors will broaden the distribution of separations via a Gaussian 
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FIG. 1: Effect of photo-z errors on the mean pairwise velocity as a function of the three-dimensional separation r (left panel) 
and on the projected mean pairwise velocity as a function of angular separation 8 (right panel). The solid line represents the 
spectroscopic sample where the positions of the SNIa host galaxies are known accurately. The dashed line corresponds to a 
scenario where the dispersion in the photo-z distribution about the true redshift is given by a z — a z o(l + z) with <r 2 o = 0.01, 
whereas the dot-dashed line represents the case when a z o = 0.02. Projected statistics vary less with cr z o, so they are less 
sensitive to systematic errors in this quantity. 



kernel with dispersion a = ^/2a z = ^/2g z q(1 + z), which gives a z — 0.056 at z = 1 for er 2 o = 0.02. To the degree 
that the assumption of small \z p2 — z p i\ is violated, the small distance error induced by this approximation remains 
negligible, as the pairwise velocity does not vary rapidly on any scales of interest. One additional caveat is that 
these relations hold only for sufficiently large angular separations, corresponding to comoving separations greater 
than approximately 5 Mpc, so that nonlinear effects due to velocities within gravitationally bound objects ("fingers 
of god") are insignificant. 

The left panel of Figure [1] shows the effect of photo-z errors on mean pairwise velocity measurements as a function 
of three-dimensional separation r. For a photo-z error <j z q — in the range 0.01 to 0.02, the overall amplitude of the 
mean pairwise velocity is suppressed by a factor of 3 to 4 at separations r < 50 Mpc/h compared to the case where all 
redshifts are known perfectly. As the separation increases, this suppression becomes less prominent. This is largely 
because the three-dimensional separations of the SNela are uncertain by an amount given by the photo-z error, which 
may be large compared to the three-dimensional distances when the separation is small. For pairs which are farther 
apart (and often have distances dominated by their transverse separations), this smearing effect has much less impact. 
The right panel of Fig. [T] shows the projected mean pairwise velocity as a function of angular separation for different 
assumed photo-z errors. Note that because of the integration along the line of sight, changing photo-z errors by a 
factor of two, from a z o — 0.01 to <j z q = 0.02, causes only a 10% to 15% change in the amplitude of the statistic. 



B. Evolution in the Luminosities of Type la Supernovae 

Evolution in intrinsic SNIa properties is one of the most important potential sources of systematic error that could 
bias estimates of cosmological parameters from pairwise velocities. For instance, the mean intrinsic luminosity of 
SNela could vary significantly over time. If this were unaccounted for, the inferred distance moduli for supernovae 
would have a systematic error whose amplitude is a function of redshift. Following Ref. [43| . we model an error in the 
distance modulus fi as 



M = Mtrue + A/X = //true + HLZ + flQZ , 



(19) 
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where fi^ and fiQ are parameters quantifying the linear and quadratic dependence of the systematic error on redshift. 
A systematic error in distance modulus propagates into an error in the inferred cosmological redshift z(p) in Eq. 
while a systematic error in the photometric redshift directly affects z mca s- 
Propagating errors through Eqs. and © gives 

Az= 0.46Ml + z) > (20) 



where we define the function 



(I 4- z) 2 

f(z) = l+ } \*> (21) 
d L (z)H{z) 

The resulting error in the line-of-sight velocity for a given supernova is 

0.46cA« 

Av ios = — - (22) 

J ^meas / 

Using the measured redshift instead of the cosmological redshift in this expression gives an error on the order of a 
few percent at redshifts of interest. 

This systematic shift can be applied directly to the estimator Eq. (TTTj) to evaluate the impact systematic errors have 
upon a given supernova velocity catalog. Alternately, we can apply this systematic error to Eq. ([S]) to get an estimate 
of the size of the resulting shift in the pairwise velocity statistic. Consider a pair of supernovae with measured redshifts 
z\ and Z2- Each of them has their three-dimensional velocity systematically shifted in the line-of-sight direction by an 
amount At>i os ; the component of this shift along the vector connecting the two galaxies is Av\ os TT t /r where n t is their 
separation along the line of sight and r is the distance between the galaxies. Their pairwise velocity gets a systematic 
shift given by 

Av(r, a) = — — — = ~p7 T - 77 — \ - 77"^ (ml + 2ziMq), (23) 

r \l + z 2 1 + 21/ r \f{z2) f{zx)J J{zi)r 

where for the last expression we have used the fact that the difference in the second expression is dominated by the 
difference in distance modulus, rather than the much smaller difference in f(z). In replacing both redshifts by z\ in 
this expression, we have assumed that the redshift difference for a given pair is small compared to unity, which will 
be the case for any pair separations for which the mean pairwise velocity is significant. For a given value of r and 
a=(l + zi) _1 , the only quantity which varies between different pairs is the line-of-sight separation term, 7rf , whose 
average over random pairs is r 2 /2. Averaging the final expression in Eq. (|23p over all pairs with a given separation 
replace ir 2 /r by r/2 and gives the systematic error in the mean pairwise velocity for pairs with comoving separation 
r and mean redshift z as 

A,(r,z)^0.23r g(z)( ^ + 2z/io) . (24) 
/ \ z ) 

For the systematic error in the projected statistic, we substitute Eq. into Eq. (fl^j) . which yields 

Av(6,z) =0.23^ <J*l+2hqz) f dn t P{n t \0, z)Je 2 d c {z) 2 + tt 2 , (25) 
J\ z ) Jo 

with P(n t \0,z) given by Eq. (|T5|) . This expression is used in Sec. I VII to estimate how small this systematic error 
must be so that it does not dominate the statistical errors in mean pairwise velocity measurements of dark energy 
parameters. 



V. STATISTICAL ERRORS 



The line-of-sight velocity for a supernova is inferred by combining a redshift measurement and a distance estimate 
obtained from a brightness measurement. Here we assume Gaussian random errors for both the redshift and brightness 
measurements, and find the resulting statistical error in the mean pairwise velocity. We also give an expression for 
the sample variance (sometimes referred to as cosmic variance) error in this quantity, which results from the fact that 
its intrinsic value in the limited volume we probe may not match the universal mean. 
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A. Apparent Magnitude and Redshift Errors 



For a given supernova, we assume normal errors of and a z on the distance modulus and the measured redshift. 
Propagating through Eq. ([5]) using Eq. (|20[) and adding the resulting errors in quadrature gives 

Sv los (z) 2 = ^2°l + / n I y v» °i ( 26 ) 

f(zy (l + zy 

In evaluating the first term, we have assumed that a z is small compared to z meas , which should be a good approximation 
for photometric redshifts of SNela [HI, [53| . This allows us to neglect the effect of errors in 2: mca s on the value of f(z). 
Note that in actual measurements the errors in photometric redshifts may be significantly non-Gaussian, requiring a 
more sophisticated treatment; here we explore Gaussian errors to give an approximate guideline for the relevant levels 
of uncertainty. 

Gravitational lensing may increase the dispersion in the measured distance moduli of SNela beyond that of intrinsic 
luminosity scatter and random measurement errors. In the weak lensing limit (convergence n <C 1), the dispersion 
due to lensing is 0, Q3, [H| 

r z W 2 (z' z) r 

a 2 cns (z) « 1.690^ jf dz' ^ ' } J dkkP{k,z'), (27) 

where W(z',z) = HodA(z')dA(z',z)/dA(z), dA{z) is the angular diameter distance to redshift z, and dA(z',z) is the 
angular diameter distance between redshifts z' and z. The quantity P{k, z) is the matter power spectrum; we evaluate 
it using the numerical fits of Smith et al. (56|. We thus have a total standard error on the distance modulus for a 
single supernova composed of three pieces: 

a l = a lbs + °SN + View ( 28 ) 

where <7 b s is the random scatter due to measurement noise and ctsn * s ^ ne intrinsic scatter in supernova intrinsic 
luminosity. Where not otherwise specified, we take ctsn — 0.1 independent of redshift, following recent estimates [43j , 
and assume <7 bs <C osn, which should be satisfied for upcoming large surveys like LSST. 

To obtain the standard error in the mean pairwise velocity, we begin by assuming that each individual line-of-sight 
velocity has a normally-distributed error with standard deviation Svi os . Then for any data bin, applying standard 
propagation of errors to Eq. (jll[) gives: 

6v cst = V25v 1o JJ2pI] > ( 29 ) 

pairs 

assuming that fractional errors in the Pij are modest; we expect this to hold, as these values can be evaluated using 
redshift distances, rather than the comparatively uncertain distance measurements that drive the uncertainty in 
individual speeds. For each pair, pfj ~ cos 2 tp, where ip is the angle between the comoving line-of-sight vector and 
the vector connecting the comoving supernova positions. This angle will be distributed randomly for each pair; the 
expected mean value of p 2 ^ over a large number of pairs is 0.5. Thus the standard error in the mean pairwise velocity 
in a particular redshift and separation bin is 

ftr*(r,a)=2-%W (30) 

where N(r,a) is the total number of pairs used to estimate the mean pairwise velocity in a given redshift bin with 
mean scale factor a and separation bin with mean separation r. 

The standard error on the projected statistic can be expressed as a sum over pairs in the same way as v(r, a), except 
the sum is over iV(#, z) pairs in a given angular separation bin about 9 instead of a given real-space separation r. The 
same calculation applies, except that now the average value of p 2 ^ for a bin in 9 will not be 0.5. For a given pair, the 
projector = TVt/r, where 7Tt is the comoving radial separation of the pair (as defined in section lHl|) . Analogous to 
Eq. (|29p , the error on the projected statistic in a bin can be written as 

/ 2\-V2 

Sv est = V2Sv lo Jj2^) • (31) 
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The sum must be evaluated by integrating over all the pairs in a given angular bin, giving 

-1/2 

(32) 

with P(ir t \9, z) given by Eq. ([TBI . 

For a bin in angle covering a range from 6\ w to #high and a mean redshift bin with a range from z\ ow to Zhigh, we 
can derive the number of pairs in this bin from Eq. ([T]). Consider a supernova at redshift z\. Any second supernova 
which lies in the angular bin will be contained in a sky region with area 27r(cos#i ow — cos#high) — ""(^high — ^low)' 
where the second expression is valid for small angles. The second supernova at redshift z 2 > z\ must satisfy z\ ow < 
{z\ + z 2 )/2 < Zhigh f° r the pair to be in the redshift bin, and c(z 2 — z 1 )/i?(z 1 ) < 7r max for the comoving line-of-sight 
separation to be less than 7r max . These conditions are equivalent to 

ftmaxH (^low) j /oo\ 

ziow -z < zi < z high and (33) 



5v(6,z) = Svi os (z) 



N(6,z) 



dir t P(ir t \6,z) 



d 2 d c (z) 2 + ir 2 



2z lo 



z\ < Z2 < mm 



2z h igh - zi, zi + 



TTmax-ff^l) 



Then neglecting the effect of any spatial clustering of supernovae, the total number in the bin is simply 



N(0\ ov!l #high; Z\ ow , Zhigh) — ^(^high — ^low) / dz 



d 2 



dz dfl 



(zi) / dz 2 



' dz dQ 



(* 2 ), 



(34) 



(35) 



where the limits on the z integrals are given in Eqs. (1331) and (I34[) : note the z 2 integral must be performed first since 
its limits depend on z\. The function d 2 n/dzdfl is just Eq. ([T]) normalized to the total number of supernovae assumed 
per unit solid angle on the sky. 



B. Sample Variance 



In addition to the measurement errors for individual galaxy velocities, there is an additional uncertainty in com- 
paring estimates of the mean pairwise velocity to models, resulting from the fact that we only sample a finite volume 
in which the realized average pairwise velocity may differ from the mean taken over the entire Universe. Here we 
give an expression for the covariance between the projected mean pairwise velocity measured in different redshift and 
angular separation bins resulting from this effect (generally referred to as sample or cosmic variance). 

Consider a mean pairwise velocity statistic binned in pair separation, r, and scale factor, a. For the three-dimensional 
mean pairwise velocity, Eq. (J7J), the sample covariance between two bins in separation and scale factor [r, a] m and 
[r, a] n for a survey volume Vq can be written as 23] 



_ 32TrH(a m )a m H(a n )a n 

^ v m 5 i &m j &n ) 



9Vh(i + e gal ( ))(i + £ gal ( 

d In D„ \ ( d In D„ 



d\na I \ d In 



dkk 2 \P(k)\ 2 j 1 {kr m )ji (kr n ). (36) 



We now integrate along the line-of-sight accounting for the photo-z errors and obtain an expression for the sample 
covariance of the projected mean pairwise velocity as a function of perpendicular separation, 

/>oo />oo 

C(9 m ,9 n ;a m ,a n ) = / dn™ / dn[ m) P(n[ m) \6 m ,a m )P(n ( ^\n[ m) ) 



o Jo 



x / d^ s / d7r t (n) P(7r t (n V„,a„)P(7r^|7r t (n) )C(r m ,r„;a m) a„), (37) 



using Eqs. $15§ and ([XT]). 

The total statistical error covariance matrix is the sum of the sample covariance, Eq. (|37|). and the statistical error, 
Eq. (EH): 

Ctotai(^m; a r«7 a n) = C{6 m , 8 n ] a m , a n ) + 5 mn 8v 2 (9 m , a m ) . (38) 

In the following Section, we use this total covariance matrix to estimate the observability of SNela peculiar velocities 
and their utility to cosmology. 
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VI. RESULTS 

A. The Signal- To-Noise Ratio of Projected Mean Pairwise Velocity Measurements 

As seen in Fig.[TJ the projected velocity statistic given by Eq. (fl"2j) is far less sensitive to photometric redshift errors 
than the non-projected pairwise velocity. We therefore will use this statistic both to estimate the signal-to-noise of 
pairwise velocity measurements and to determine the resulting constraints on cosmological parameters. The simple 
pairwise velocity should yield comparable or better constraints in the limit of small photometric redshift errors, but 
the results will be more sensitive to a z . 

Figure [5] shows the signal-to-noise ratio per angular bin for measurements of the projected mean pairwise velocity 
as a function of angular separation 9, for our fiducial survey giving 3 x 10 5 total host galaxies over 300 square degrees 
of sky, and a distance modulus scatter for each host galaxy of ctsn =0.1 plus the scatter due to lensing magnification. 
Pairs are binned in 6 redshift bins equally spaced between z = and z = 1.2. For each redshift bin, 10 bins in angle 
are used, equally spaced for angles ranging from = up to the angle subtended by our maximum pair separation of 
100 Mpc at the mean redshift for the redshift bin. The maximum angle considered therefore decreases as the redshift 
increases, causing the curves in Fig. [5] to truncate at differing values of 9. 

The mean pairwise velocity is detectable at a wide range of angular separations and redshifts. The top left panel of 
Figure [5] shows that such a measurement with a photometric redshift error of a{z) = 0.01(l + z) yields a signal-to-noise 
ratio between 2 and 9 over a range in angular scales for all but the most extreme redshift bins with z > 0.8. The 
redshift distribution of observed SNela peaks around z — 0.5 in our LSST-likc model, so the closer we get to that 
redshift range, the more host galaxy pairs we average over and the better we can measure velocity statistics. Note 
also that although the number of pairs increases at larger separation, the amplitude of the mean pairwise velocity 
decreases, yielding an overall decrease in the signal-to-noise for bins with larger separations. The top right panel of 
Fig. [5] shows the signal-to-noise ratio for a photometric redshift error of a z = 0.02(1 + z). After this doubling of the 
photometric redshift error, the signal-to-noise decreases by around 30%. Even for a z = 0.03(1 + z) (lower left panel 
of Fig. [2]), we still reach a signal-to- noise of around 3 for the redshift bins at z — 0.5 and z — 0.7. 

The lower right panel shows a best-case scenario, assuming that spectroscopic redshifts are obtained for each 
supernova host galaxy; for simplicity, we define a spectroscopic redshift to have a z = 0.001(1 + z). This redshift 
error is generally obtainable only from spectroscopy of the hosts (rather than the SNe themselves), primarily because 
of the large breadth of SNIa spectral features, but also due to the peculiar velocities of SNe with respect to their 
galaxy's center, which can reach a few hundred km s _1 . Spectroscopic redshifts for large samples of hosts (though 
likely not all, since many will be fainter than the SNe) would be quite feasible with a 5000-fiber, large field of view 
multi-object spectrograph like that currently proposed for the BigBOSS project 57]. If supernova samples cover 300 
square degrees, as assumed above, a minimum of 43 BigBOSS pointings would be required to cover this sky region, 
yielding more than 200,000 redshifts; larger samples can be obtained by revisiting each pointing with different fiber 
placements. The proposed BigBOSS survey would use the Kitt Peak 4-meter telescope for only 100 nights per year; 
such a supernova project would require only a small fraction of the remaining time available. In this "spectroscopic 
limit," the signal-to-noise in measuring the mean pairwise velocity generally improves by around a factor of two 
compared to the o z = 0.01(1 + z) case. 

B. Parameter Space and Formalism 

Now we investigate the constraints on dark energy parameters from a SNIa projected mean pairwise velocity 
measurement, and assess the complementarity of these constraints to performing the luminosity distance test based 
on the same data. For the sake of simplicity, we perform a Fisher matrix analysis similar to those in Refs. fl7l. [23j| . In 
order to compute constraints on f?A, wq and w a , we marginalize over the remainder of the parameter space, consisting 
of the parameters Af, ns, and h. We also treat a Z Q, describing the photometric redshift dispersion, as a parameter 
since the binned mean pairwise velocity signal depends on this quantity. 

In addition to the marginalized constraints on wq, and w a , we quantify the additional constraining pow er of 
pairwise velocities by evaluating the quantity [a(w p )a(w a )]~ 1 for comparison to the DETF summary tables [43|]. We 
refer to this as the "Figure of Merit" (FoM) for convenience, although in the DETF report this term refers to a slightly 
different quantity (the inverse area of the 95% confidence limit ellipse in the w p — w a plane) which is proportional 
to [a(wp)a(w a )]~ 1 . The derived parameter w p is the equation of state at the "pivot" (i.e. best-constrained) redshift, 
defined as w p = w + (1 - a p )w a with a p = 1 + [F~ 1 } WoWa /[F~ 1 ] WaWa . 
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The Fisher matrix for the projected mean pairwise velocity can be written as 



dv(m) 



C toU mn ) 



dv(n) 



(39) 



dp, 



dpp 



where we have abbreviated the projected mean pairwise velocity in the nth angular separation and redshift bin as 
v(n), Ctotai(wT') is the total covariance matrix between bins m and n given by Eq. (|38p , and p a indexes the parameters 
in the vector p. The Fisher matrix provides a local estimate of the parameter covariance, so the standard error on 
parameter p a marginalized over the other parameters is u{p a ) = [-F _1 ] Qa (no summation implied). 

Prior constraints on any of the parameters p which are normally distributed are simple to incorporate. If parameter 
p a has a Gaussian prior with standard error er Q , we simply add the diagonal matrix diag(l/cr^,) to the Fisher matrix 
F a p. Priors with non- normal statistical distributions require a more detailed statistical framework rather than a 
simple Fisher matrix approximation. 



In computing constraints on the dark energy parameters wq, and w a , we first assume a reasonable calibration 
spectroscopic sample of 1500 SNela, comprising 250 supernovae in each redshift bin spread uniformly over the 6 
redshift bins spanning < z < 1.2. The fractional error on the photo-z dispersion, 5o~ z o/o~ z o in this case is around 
1/V500, or approximately 5% (assuming Gaussian errors). We therefore incorporate a Gaussian prior on o~ z o centered 
on the true value and with a = 0.05ct2o; however, as we show below in Fig. 31 the pairwise velocity statistic is relatively 
insensitive to the choice of a prior on <t z q, so this choice should not significantly affect our results, even if the actual 
error on o~ z q is much larger. 

We compute the standard errors obtainable on the dark energy parameters using a range of supernova distance 
modulus dispersions ctsn and photometric redshift dispersions <t z q. We consider three possible values of the intrinsic 
supernova absolute magnitude dispersion given by ctsn = 0.05, osn = 0.1, and ctsn = 0.2. For each value of ctsn, we 
explore four possible values of photometric redshift dispersion, (t z q = 0.001 (the "spectroscopic limit"), a z o = 0.01, 
(j z o = 0.02, and <r z o — 0.03. The optimistic but reasonable supernova luminosity distance test assumed in the DETF 
report corresponds to ctsn = 0.1 and ct z o = 0.01 so these choices constitute a sensible baseline for comparison to other 
techniques. 

The strength of the dark energy constraints obtained is relatively sensitive to the amount of prior information 
assumed. First, we can make the same assumptions used by the Dark Energy Task Force (43[. They assume constraints 
on all parameters (including covariances) at the level expected for measurements of the microwave background power 
spectrum by the Planck satellite. For this, we employ the Planck Fisher Matrix provided by the DETF. In addition, 
DETF assume a 11% Gaussian prior on the value of h [Hj]. Note that a spatially flat universe is not assumed. We 
also assume no systematic error on either redshift or distance modulus measurements; limits on these systematics 
required to attain the statistical error levels presented here are discussed below. The results are given in Tabic [I] 

For the nominal DETF survey case, mean pairwise velocities give a standard error on wq of ct(wo) = 0.45 and 
a standard error on w a of o-(w a ) = 0.98. This constraint on wo is comparable to the DETF Stage-IV constraints 
from ground-based optical baryon acoustic oscillations or galaxy cluster counts, while not as good as those from 
Stage-IV supernova luminosity distances. For w a , mean pairwise velocity constraints are significantly better than the 
optical survey-based BAO projection; slightly weaker than the pessimistic BAO projections for space-based or radio 
observations and for the optimistic galaxy cluster projection; and halfway between the optimistic and pessimistic 
DETF supernova luminosity distance projection. However, all of these methods trail the Stage IV weak lensing 
projections in constraining power. 

Among the dark energy probes resulting from a large ground-based optical survey like LSST, mean pairwise velocities 
compare well with both the baryon acoustic oscillation and the supernova luminosity distance probes [13] ■ To quantify 
this, we consider the improvement in dark energy parameters obtained by adding the mean pairwise velocity probe to 
the supernova luminosity distance probe resulting from the same sample. The mean pairwise velocity can be measured 
using the supernova data from a large survey telescope with little additional cost compared to simply constraining 
dark energy using the resulting supernova Hubble diagram. 

Figure [3] shows joint constraints on the dark energy parameters combining projected peculiar velocity measurements 
and the SNIa luminosity distance test, using the same priors as TableQ] The left panel shows the lex constraint in the 
wq — J1a plane and the right panel shows the constraint in the w a — plane, after marginalizing over the remainder 
of parameter space. Incorporating peculiar velocity information significantly reduces the size of the ellipses in the 
dark energy parameter space: the marginalized constraint on improves by a factor of 1.7, on wq by a factor of 
1.2, and on w a by a factor of 1.5, giving an overall improvement in the Figure of Merit by a factor of 1.8. (As a 



C. Statistical Constraints on Dark Energy Parameters 
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O"z0 


= 0.001 




O"z0 


= 0.01 




(TzO 


= 0.02 




<7>0 


= 0.03 




C"SN 


o-(Oa) 


a (w ) 


cr(w a ) 


o-(Oa) 


cr(u)o) 


cr(wa) 


<r(n A ) 


cr(wo) 






cr(wo) 


<x(w a ) 


0.05 


0.016 


0.22 


0.36 


0.032 


0.28 


0.48 


0.05 


0.43 


0.89 


0.059 


0.78 


1.7 


0.1 


0.039 


0.30 


0.55 


0.056 


0.45 


0.98 


0.094 


0.62 


2.52 


0.13 


1.46 


2.81 


0.2 


0.051 


0.42 


0.72 


0.074 


0.59 


1.84 


0.18 


1.36 


4.71 


0.24 


2.87 


6.14 



TABLE I: Dark energy parameter constraints derived from mean pairwise velocity statistics. Photometric redshifts are assumed 
to be normally distributed about the true z, with <r z — <r z o(l + z). We show results for <r z o = 0.001, 0.01, 0.02, and 0.03, and 
three different values for the uncertainty in supernova distance moduli, <jsn = 0.05, 0.1 and 0.2. The fiducial values of the 
dark energy parameters are Qa = 0.75, Wo = —1 and w a = 0. We assume zero systematic errors related to SNIa evolution, i.e. 
Mi = — 0. We assume the same priors used in the DETF report: Planck satellite priors from its projected measurement of 
the microwave background power spectrum (using the Fisher matrix supplied by the DETF) and a Gaussian prior on h with a 
standard error of 11%. This table does not assume a fiat spatial geometry for the universe. 

point of comparison, corresponding constraints with no priors from other measurements are included in Fig. 2J) Note 
that unlike the case of peculiar velocity measurements, the constraints derived from the SNIa luminosity distance 
are sensitive to the error in mean redshift of a bin and hence the cosmological constraints derivable from the SNIa 
luminosity distance depend much more on the amount of prior knowledge of the photo-z distribution [TtI |59| , as well 
as being much more sensitive to intrinsic SNIa luminosity evolution. 

We have also considered statistical dark energy constraints from a more constraining, but still realistic, set of 
priors. In particular, a measurement of the Hubble parameter based on an improved, NGC 4258-calibrated distance 
ladder with an estimated overall error of 5% has recently been reported [60[. Furthermore, requiring the dark energy 
probe itself in combination with microwave background data to determine the geometry of the universe is likely 
overly restrictive. Measurements of the baryon acoustic oscillation scale from the Sloan Digital Sky Survey Data 
Release 7, combined with WMAP 5-year data, give a constraint on the curvature parameter of Slfc = —0.013 ± 0.007, 
even for a very general cosmological model which allows both a nonfiat universe and a value of wq different from -1 
(6l| . Additionally, since a flat universe is an unstable fixed point for standard cosmological evolution, we have an 
overwhelming theoretical prejudice for fife = to high precision. Therefore, a prior assumption of a flat universe is 
both reasonable and strongly suggested by data. 

Table |TT] gives the standard errors on the dark energy parameters for a flat universe, with gaussian priors for h (5%), 
Af (5%), and n$ (1%), the latter two being current limits from WMAP 7-year data (62[. The assumption of a flat 
universe and a tighter prior on h lead to much stronger dark energy constraints than do DETF priors. With these 
priors, a DETF-assumcd supernova sample with ctsn = 0.1 and a z o = 0.01 gives a measurement of f^A, wq and w a 
with standard errors of 0.024, 0.27, and 0.41, respectively, using the mean pairwise velocity alone. For comparison, 
the constraints on wq for Stage IV experiments computed in the DETF report (but for the original set of priors) are 
worse for clusters, comparable for baryon acoustic oscillations, and better for the supernova Hubble diagram. Our 
constraint on w a , on the other hand, is better than for any of the Stage IV experiments aside from the optimistic 
weak lensing scenarios. 

Of course, the constraining power of other probes will also increase with the more restrictive set of priors we assume 
in Table HU This makes a direct comparison with these other methods beyond the scope of this paper. Our primary 
point is that under optimistic, but reasonable, assumptions, SNIa peculiar velocities can be useful by themselves and 
at the very least can serve as a valuable complementary probe and cross-check for systematic errors, while requiring 
little additional investment. However, note that Tabic [TTI shows that broader photo-z distributions and/or larger 
intrinsic SNIa dispersions can quickly diminish the returns on SNIa peculiar velocities. 

This calculation also suggests the potential constraining power of pairwise velocity statistics from future survey 
observations. If a large photometric supernova survey were combined with follow-up spectroscopic redshifts for 
supernova host galaxies, the standard error in the redshift could be reduced by a factor of 10 to a zQ — 0.001, 
corresponding to the first column of Table HTl In this case, the error on Wo shrinks to <r(wo) ~ 0.10 and the error 
on w a is nearly a(w a ) ~ 0.16. Understanding Type-la supernovae well enough to push <tsn down by a factor of 
2 to ctsn = 0.05 would reducethe error on w a by another factor of two, to 0.08. Few other proposed probes have 
comparable potential to constrain uu a . 
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a z0 = 0.001 


cr z0 = 0.01 


<Tzo = 0.02 


CTzo = 0.03 


C"SN 


ct(Oa) ct(wq) a(w a ) 


ct(Oa) a(wo) o(w a ) 


ct(Oa) cr(iuo) o~(w a ) 


<t(Qa) <t(wo) &(w a ) 


0.05 
0.1 
0.2 


0.004 0.047 0.08 
0.009 0.10 0.16 
0.022 0.28 0.41 


0.012 0.16 0.23 
0.024 0.27 0.41 
0.061 0.63 0.89 


0.024 0.34 0.49 
0.046 0.62 0.94 
0.1 1.48 1.93 


0.049 0.65 1.03 
0.1 1.42 1.96 
0.2 2.86 4.18 



TABLE II: Same as Table H but we assume = and instead of Planck priors we assume Gaussian priors with standard 
errors of 5% on h and Af and 1% on ns (comparable to errors from current measurements). 



P-L = PQ 


Q.A WO W a 


0.01/V2 
0.03/^2 
0.05/\/2 


10.2% 4.1% 20.2% 
20.3% 10.0% 41.0% 
40.9% 37.5% 80.1% 



TABLE III: The ratio of the parameter bias due to systematic error to the statistical uncertainties on these parameters. A 
photometric redshift distribution with dispersion a z = 0.01(1 + z) is assumed. We assume ph = pq to compute the systematic 
bias. Then we set = PQ = and compute the statistical uncertainty and report the ratio of systematic bias to statistical 
errors Ap/a p , where p — Oa, Wq, or w a ■ We assume fife = and assume Gaussian priors with standard error of 5% on h and 
Af and 1% on n s . 

D. Systematic Error in Distance Modulus 

The potential statistical sensitivity of any dark energy probe can only be realized if systematic errors can be 
controlled to a level where their effect on cosmological parameters is small compared to the statistical errors. For 
the supernova data set considered here, systematic errors may effect both observables: the distance modulus and the 
photometric redshift. This section considers distance modulus systematics, while the following section analyzes the 
effect of redshift errors. 

Section TlV Bl gives a simple phenomenological model for the effect of SNIa evolution with redshift, in terms of the 
parameters ph and pq. The resulting systematic error on cosmological parameters induced by this systematic error 
can be estimated using a Fisher matrix approach. The bias in parameter p a can be written as 

Sp a = J2i F ^U E A ^)C^Umn)^ (40) 

j3 m,n P 

where Av, obtained by substituting Eq. for v(r,a) in Eq. (fT2)l . is the systematic shift in the observable v due to 
the systematic error characterized by nonzero values of p L and fig . 

We calculate the bias in each parameter due to SNIa evolution assuming a photometric redshift distribution with 
spread a z — 0.01(1 + z) and the evolution model given by Eq. (fTT)|). We can then compare the systematic bias with the 
statistical errors on dark energy parameters assuming pl — PQ — 0, as computed in Table [TT1 The ratios of the bias of 
the dark energy parameters to their statistical errors are reported in Table III for several representative choices of /i L 
and [iq. For reference, DETF took evolution in SNIa luminosity with hl = Pq = O.Ol/v 7 ^ as their optimistic scenario. 
We find that the maximum bias incurred in 51a and Wq is less than 40% as large as the statistical error on these 
parameters as long as /il = Mq < 0.05/V2 (five times larger than the DETF optimistic systematic error). For w a , the 
systematic bias is 40% of the statistical error for [i L = fj,Q < 0.03/^/2, and increases to 80% of the statistical error for 
Ml = Pq = 0.05/V2. If the actual unrecognized evolution of SNIa luminosity is similar to that assumed in the DETF 
report, the resulting systematic bias in dark energy parameters should be insignificant compared to the statistical 
error. Note that these comparisons conservatively use the statistical error incorporating our more restrictive prior 
than in the DETF report. The larger statistical errors with the DETF priors admit substantially larger systematic 
errors. 



E. Systematic Errors in Photometric Redshifts 

We have also tested how a possible bias Az v in the photo-z distribution might impact the dark energy constraints 
obtainable from pairwise velocity statistics. If Az p is not a strong function of redshift (i.e., it does not vary considerably 
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<j{<Tz0) 


ct(Oa) a(wo) o(w a ) 


no prior 
prior (100% error in cr z o ) 
prior (zero error in cr z o ) 


0.046 0.29 0.62 
0.023 0.25 0.37 
0.018 0.21 0.31 



TABLE IV: Impact of prior information about photo-z distributions on the dark energy parameter constraints derived from 
mean pairwise velocity statistics (calculated in the absence of complementary cosmological probes). The photometric redshift 
of an SN is assumed to be normally distributed about its true value with o z = a z o(l + z); for our standard scenario we take 
tJzO ~ 0.01. The fiducial values cosmology considered has S1a = 0.75, wo = —1 and w a = 0. We assume Qk = and assume 
Gaussian priors with standard error of 5% on h and and 1% on n s . 

within one of our redshift bins with width 5z ~ 0.2), then the bias affects both galaxies in each pair in approximately 
the same manner. The mean pairwise velocity relies on the difference between the two velocities so nearly all of the 
effects of a photo-z bias tend to cancel. The residual is a small misestimation of the location of the redshift bin, which 
translates into a small error in cosmological parameters. For example, assuming a bias in photometric redshifts of 
Az p ~ 0.002(1 + z) degrades the constraints on cosmological parameters by less than 2% of the statistical errors. This 
stands in stark contrast to the strong dependence of the luminosity distance test on photometric redshift biases (e.g., 
^,[5{|) and the similar sensitivity of probes such as weak gravitational lensing to biased photometric redshifts (e.g., 

The signal we measure, the redshift-binned projected mean pairwise velocity Eq. (|18[) . depends on the scatter in 
photometric redshifts so we also must estimate the systematic error due to uncertainty in the photometric redshift 
dispersion. We assume that the distribution of the difference between photo-z's and spectroscopic redshifts is a 
standard normal distribution; in reality this distribution is likely more complex. The results here are a simple 
effective model for the distribution of photometric redshifts. 

Figure [4] and Table IIVI show marginalized statistical constraints on dark energy parameters from mean pairwise 
velocity only, under three strongly different assumptions regarding the photometric redshift error. The blue (gray) 
and the black shaded regions show the two extreme cases. The blue shaded area shows the la constraint when we 
assume no prior knowledge of the uncertainty in the photo-z error and allow a Z Q to be determined from the same data 
used to constrain cosmology. The black region indicates the constraints when ct z q is known exactly. We emphasize 
that this does not mean that the photometric redshift is equal to the true redshift. There is still a non-negligible 
dispersion in photometric redshifts in this case; however, we have assumed that the photometric redshift distribution 
is well understood, perhaps due to calibration with several thousand spectra [l7|. The red (light shaded) region 
represents the case when the prior on a z o is a Gaussian centered at the true value with sigma equal to its fiducial 
value, <7 z q = 0.01. 

Constraints on wo, w a and fl\ change by only about 10% between the case where <r z o is uncertain at the 100% level 
and one where we assume a perfectly-calibrated photometric redshift distribution. This results from the fact that the 
mean pairwise velocity is proportional to the redshift difference between galaxies in a pair, but photometric redshift 
errors do not correlate with the velocity we are trying to measure. Fig. 3] shows that even weak prior knowledge of 
the photo-z distribution yields constraints comparable to a scenario where the photo-z error distribution is known 
exactly. 

VII. DISCUSSION AND FUTURE PROSPECTS 

With vastly increased numbers of Type la supernova detections on the horizon, a new statistical probe of dark 
energy using supernova peculiar velocities will be possible. The Dark Energy Task Force, when considering future 
supernova measurements, made the optimistic but reasonable assumptions that individual supernovae will have a 
photometric redshift determined with a standard error 0.01 at redshift z = 0, and a distance modulus determined 
with an error of 0.1. If these levels are attained for the nominal 3 x 10 5 SNela which will be detected in a targeted 
supernova survey area by the LSST, the resulting dark energy constraints from the mean pairwise velocities of these 
supernovae are interestingly good, comparable to projections for a variety of Stage IV techniques. In particular, 
pairwise peculiar velocities alone give a slightly stronger dark energy constraint as the optimistic projection for an 
optical baryon acoustic oscillation probe, and constraints which are towards the optimistic ends of the galaxy cluster 
abundance and optical supernova Hubble diagram probes. 

Having another independent method for constraining dark energy is invaluable, since all of these measurements 
will likely be limited by systematic error control. Comparison of inferred dark energy parameters from multiple 
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independent experiments is even more important than the combined statistical power of multiple measurements. In 
addition, the science return from mean pairwise velocity measurements comes essentially "for free," as it uses the 
same data sets from which supernova luminosity distance measurements will be built. 

An extension of these observations which can significantly improve the strength of dark energy constraints is the 
addition of spectroscopic redshifts. Surveys like Pan-STARRS and LSST will provide only photometric redshifts, and 
the sheer number of objects they observe makes obtaining spectroscopic redshifts for even small subsets of the total 
objects a massive challenge. Given the numbers used in this paper, LSST will detect on the order of 100 new SNela per 
night of the survey. The supernovae themselves are transient and widely spread over the survey area, limiting the total 
number which can be observed simultaneously. Obtaining immediate redshift follow-up for all of these objects would 
be a large logistical challenge, even if a dedicated telescope were available. An alternative is to obtain redshifts of host 
galaxies after their SNe have faded; this can be done much more efficiently, as many host galaxies in a particular field 
of view could be targeted simultaneously with multi-object spectrographs. Baryon acoustic oscillation observations, 
in particular, are pioneering the development of very large fiber spectrographs which can obtain thousands of redshifts 
simultaneously. The challenge for this strategy is that many of the host galaxies will be at redshifts above 0.5, and 
many of the host galaxies themselves are dim compared to their supernovae. Host galaxy followup would require a 
large telescope and large amounts of observing time. The spectrograph proposed for the BigBOSS survey, which has 
been designed to obtain high-throughput spectroscopy of 5000 galaxies at a time over a 7 square degree field of view 
using the 4m Mayall telescope at Kitt Peak and the Blanco telescope at CTIO [I?! , would be well suited for this task. 
Such a spectroscopic survey of supernova hosts would be a major undertaking, but could lead to highly competitive 
dark energy constraints and leverages instruments and data already planned for other purposes. 

Another possible avenue for improvement is better standardization of SNIa intrinsic luminosities. Here our baseline 
assumption, along with the DETF, is that SNIa distance moduli will be known with a standard error of around 0.1. 
It is an open question whether we eventually will understand the SNIa explosion mechanism in enough detail, and 
have sufficient observational information, to model some portion of this scatter and reduce the effective random error. 
Magnification due to gravitational lensing provides an additional source of scatter, which can be partly understood due 
to its strongly non-Gaussian distribution, but for our nominal model survey we are not limited by lensing scatter. The 
marginalized constraint on w a , the most challenging dark energy parameter to measure, can be improved by a factor 
of two if the scatter in the intrinsic supernova distance modulus is halved. A mean pairwise velocity measurement for 
3 x 10 5 supernovae with spectroscopic redshifts and an intrinsic distance modulus scatter of 0.05 would constrain w a 
with a standard error of 0.08 using our set of current prior constraints. 

The statistical power of any given dark energy measurement is only half of the story, as all of these measurements 
are likely to be heavily dependent on systematic error control. Because of its nature as a differential measurement, 
the mean pairwise velocity technique offers favorable prospects for controlling systematic errors. Differential measure- 
ments have long been exploited in measurements of the cosmic microwave background fluctuations precisely for their 
systematic error advantages. In particular, we have demonstrated that several obvious systematic error sources are not 
likely to dominate the dark energy constraints. First, uncertainty about the level of scatter in photometric redshifts 
about their true values has only a weak effect on dark energy constraints, and mild priors obtainable from modest 
spectroscopic calibration efforts give results that are nearly the same as exact knowledge of the photometric redshift 
scatter. We have not considered non- Gaussian errors in photometric redshifts, but any scatter which is characterized 
at the levels of the normal errors considered here is unlikely to induce any significantly larger systematic errors. 

Second, a bias in the photometric redshift distribution has very little effect on our constraints, as long as the bias 
varies slowly with redshift, because a constant redshift bias doesn't affect the pairwise velocities. This is in marked 
contrast to both the supernova luminosity distance and weak lensing techniques. Both of these widely discussed 
routes to dark energy constraints are very sensitive to photometric redshift biases (l7l l59| , where a redshift bias can 
mimic a shift in dark energy parameters. Third, a systematic error in distance modulus due to unrecognized evolution 
in mean supernova luminosity with redshift will be a small effect provided the magnitude of the error is within a 
factor of 3 of that considered in the DETF report. This potential source of error can also be addressed by testing the 
rich information in supernova spectra and time series at different redshifts for any evidence of evolution in intrinsic 
supernova properties. While detailed modeling of potential systematic errors is required to understand any particular 
experiment, it is plausible that the systematic errors associated with mean pairwise velocities will be substantially 
less severe than other leading techniques for probing dark energy. 

We also note that mean pairwise velocities can be used to constrain gravitational explanations for the accelerating 
expansion of the universe. This technique has the advantage of probing structure growth over a wide range in redshift, 
while also being sensitive to the expansion rate; the comparison between these two quantities is the key to constraining 
alternate gravity models [64-66]. Pairwise velocities from a much smaller sample of galaxy clusters, with more precise 
velocities obtained via the kinematic Sunyaev-Zeldovich effect, have already been shown to offer potentially interesting 
constraints on modifications of gravity [l4[ . 

The pairwise velocity statistic offers a particularly simple route to a probe of modified gravity. In linear perturbation 
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theory, the evolution of the growth factor D(a) is given to a very good approximation by d\nD/d\x\a = fi (a) 7 .where 
7 is nearly a constant and takes the value 7 ~ 0.55 for general relativity [3]; see Ref. Juszkiewicz et al. [67j for 
a highly accurate approximation to D(a). Other gravitation theories can have different values of 7; for example, 
DGP gravity [68| has 7 = 0.68 [8]. Examining Eq. ([7]), we see that the mean pairwise velocity (on the left side) 
depends linearly on dlnD/dlna, as well as H(a), the (linear regime) galaxy bias factor, and correlation function 
information. Other cosmological tests, such as the supernova Hubble diagram, will directly constrain H(a), while 
correlation functions will be measurable directly from the data set used. The linear clustering bias of host galaxies 
can be constrained in a number of ways; e.g., by direct comparison of galaxy correlation functions to the matter power 
spectrum derived from gravitational lensing; with galaxy three-point correlation function s I6QB o r angular bispectra 
[70( 1 ; or (if a large spectroscopic sample is available) by combining redshift-space distortions |19h21| with mean pairwise 
velocity statistics. Assuming these other quantities will be measured with errors which are small compared to our 
velocity errors, a measurement of v (r, a) will provide an estimate of d In D/d In a in several bins in a, which can then 
be used to constrain 7. 

If for each bin in a we have independent measurements of v(r, a) in five radial bins with a signal-to-noise ratio of 
around 3 in each bin (see Fig. [2]), then the amplitude of the function v(r, a) can be constrained with a fractional error 
of around 0.33/-\/5 or 0.16. Assuming this is the dominant error, the fractional error on din D/dlna is also around 
0.16. By propagation of errors, the resulting error on 7 is then 0.16/ lnf2 m (a), and hence ranges from 0.15 to 0.45, 
depending on the redshift bin. This would give approximately a 25% to 80% measurement of 7 in each redshift bin 
(assuming 7 takes its general relativistic value), providing a significant constraint on many theoretical alternatives 
to general relativity. With spectroscopic redshifts, these constraints would improve by a factor of three due to the 
increase in signal-noise ratio in each angular bin; then the best redshift bin alone might provide a 10% measurement 
of 7, comparable to projected constraints from weak lensing [l2j ]. Prospects for constraining modified gravity with a 
large supernova survey will be explored in more detail elsewhere. 

Dark energy is simultaneously one of the most important problems in physics today, and one of the most elusive 
to address observationally. Mean pairwise velocities extracted from a large survey of SNela can provide an important 
arrow in the dark energy quiver and should be considered alongside any of the other methods now being actively 
pursued. If simply piggybacked on existing plans for supernova luminosity distance tests, pairwise velocities offer 
independent dark energy constraints which are competitive with other methods. If augmented by spectroscopic 
redshift followup observations, pairwise velocities alone may provide important constraints on dark energy, with 
constraints on w a of 0.1 or better. Perhaps most importantly, this technique provides not only statistical power but 
potentially strong control of systematic errors. Additionally, it allows tests of the nature of gravity which cannot be 
obtained using distance measurements alone. We anticipate that Type la supernova peculiar velocity statistics will 
be in the vanguard of dark energy constraints over the coming years. 
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FIG. 2: The signal-to-noise per angular bin for the projected mean pairwise velocity for different redshift bins as a function of 
angular separation 8, for a catalog with 3 x 10 5 total host galaxies over 300 square degrees of sky, and a distance modulus scatter 
for each host galaxy of ctsn =0.1 plus the scatter due to lensing magnification. The pairs are binned by their photometric 
redshifts and the central values of the redshift bins are shown; each redshift bin has a width of Az = 0.2. The maximum angle 
for each redshift bin corresponds to the angle subtended by 100 Mpc at the bin's mean redshift; the range in angles from to 
the maximum angle is divided into 10 angular bins. The top panels show the signal-to-noise for a photo-z normal error given 
by Oz = o z o(l + z), where a z o = 0.01 (left) and a z o = 0.02 (right). The lower left (right) panel shows the signal-to-noise for 
a z o = 0.03 (left) and the "spectroscopic" limit with a z o = 0.001 (right). 
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FIG. 3: Joint dark energy constraints obtainable from an LSST-like future supernova survey, combining constraints from the 
supernovae luminosity distance test, priors on Hubble parameter [58j and constraints from Planck with those obtainable from 
supernova velocity statistics. A flat universe is not assumed. A fiducial value for the photometric redshift error of a z o = 0.01 
is assumed, with a Gaussian prior on a Z Q of 5%. The Fisher matrices for the Planck priors and the SNIa luminosity distance 
priors are obtained from the DETF report [13] . The luminosity distance Fisher matrix represents the DETF LSST supernovae 
optimistic (LST-o) survey. The blue (dark) shaded region shows the constraint obtainable from the supernova luminosity 
distance test only. The red (grey) region shows the constraint when priors from a Planck survey are combined with the distance 
test. The green (innermost, light shaded) region shows the joint constraint combining the distance test, a Planck prior and 
mean peculiar velocity measurements for the SNIa host galaxies. 
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FIG. 4: Dark energy parameter constraints from mean pairwise SNIa velocities (in the absence of complementary cosmological 
probes), for the same survey as in Fig. [2] The left and the right panels show the la contour in the u>o — Qa and the w a — 
planes. Photometric redshift errors of <r z = 0.01(1 + z) are assumed. The black shaded region represents the case when 
the photo-z distribution is known accurately. The red (light shaded) region applies a prior such that the uncertainty in the 
photometric redshift error, 8a z o, is equal to the value of <r z o; this is highly conservative. The blue (gray) shaded region shows 
the constraint when we have no prior knowledge about the photo-z error distribution (e.g., zero supernovae with spectroscopic 
redshifts) . 



